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In this paper we summarize some aspects of the wind accretion theory in high mass X-ray binaries 
hosting a magnetic neutron star and a supergiant companion. In particular, we concentrate on the 
different types of interaction between the inflowing wind matter and the neutron star magneto- 
sphere that are relevant when accretion of matter onto the neutron star surface is largely inhibited; 
these include inhibition by the centrifugal and magnetic barriers. We show that very large lumi- 
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transients (SFXTs), a new class of high mass X-ray binaries in our galaxy recently discovered 
with INTEGRAL. According to this interpretation we argue that SFXTs which display very large 
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detect and study accreting magnetars in binary systems. 
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1. Introduction 

Supergiant Fast X-ray transients (SFXTs) are a new class of high mass X-ray binaries (HMXB), 
recently discovered with INTEGRAL. These sources are observed to exhibit sporadic outbursts, 
lasting from minutes to hours, with peak X-ray luminosities between ~10^^ and 10^^ erg s^^ 
(Sguera et al., 2005). No firm orbital period measurement has been obtained yet'. A recent list 
of confirmed (~5) and candidate (~6) SFXTs is given by Walter & Zurita Heras (2007). Be- 
tween outbursts, SFXTs remain in quiescence with luminosities in the range ~10^'-10''^ erg s^^ 
In some cases, very high peak-to-quiescence X-ray luminosity swings (factor of ~10^-10^) were 
seen on timescales comparable to the outburst duration. Some SFXTs showed also flare-like ac- 
tivity at intermediate luminosity levels (e.g., In't Zand, 2005). Optical identifications of SFXTs 
show that these sources are associated to OB supergiant companion stars (see e.g., Walter & Zu- 
rita Heras, 2007, and reference therein). These stars have typically a mass of M*~3OM0, opti- 
cal luminosity of log (L*/L0)~5-6, mass loss rate of Mw=10^^-10^^ Mq yr^\ wind velocity of 
Vw~ 1000-2000 km s^^ and are persistent soft X-ray sources with luminosity around ~10^^ erg s~' 
(Cassinelli et al., 1981). 

It is widely believed that SFXTs contain neutron stars (NS) sporadically accreting matter from 
a supergiant companion, and in the prototypical SFXTs XTE J1739-302 (Sguera et al., 2006) and 
IGR J16479-4514 (Walter et al., 2006) some evidence has been reported for spin periods in the 
~ 1000-2000 s range. Proposed models for SFXTs generally involve accretion onto a NS immersed 
in the clumpy wind of its supergiant companion (In't Zand, 2005; Walter & Zurita Heras, 2007). 
However, if accretion onto the collapsed object of SFXTs takes place both in quiescence and out- 
burst, then the corresponding X-ray luminosity swing, typically a factor of ~10'*-10^, would require 
wind inhomogeneities with a very large density and/or velocity contrast (according to the standard 

'Only IGR Jl 1215-5952 showed recurrent flaring activity, with a periodicity of ~165 d. This is interpreted as 
outbursts from a systems with an unusually long orbital period (Sidoli et al., 2007). Thus Walter & Zurita Heras (2007) 
excluded these sources from their SFXT list. 
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Figure 1 : Schematic view of a magnetized NS interacting witli tlie inflowing matter from its supergiant companion. All the regimes 
described in the text are shown, together with the relative position of the magnetospheric radius (solid line), the corotation radius (dashed 
line), and the accretion radius (dotted line). A wavy solid Hne is used when the magnetospheric boundary at the magnetospheiic radius 
is Kelvin-Helmholtz unstable. In the supersonic and subsonic propeller regime convective motions at the base of the atmosphere are 
represented with small eddies. 



wind accretion, the mass capture rate onto the NS scales like M^v"^ Davidson & Ostriker, 1979). 
Models involving accretion of extremely dense clumps are still being actively pursued (Negueruela 
et al., 2008). The requirement on the density and/or velocity contrasts in the wind can be eased 
if there is a barrier that remains closed during quiescence, halting most of the accretion flow, and 
opens up in outbursts, leading to direct accretion. After reviewing briefly the theory of wind accre- 
tion in HMXBs (§ ^ and ^), we apply barrier scenarios to SFXTs in § ^ and conclude that if SFXTs 
host slowly rotating NSs (spin periods of several hundreds to thousands seconds), then they must 
possess magnetar-like fields (^lO^'^-lO^^ G). 



2. Stellar wind accretion 



We investigate here the conditions under which a magnetized NS can accrete matter from the 
wind of a massive companion. In the theory of wind accretion in HMXBs, the following radii are 
defined (e.g., Iflarionov & Sunyaev, 1975; Stefla et al., 1986): 

- The accretion radius, Ra is the distance at which the inflowing matter is gravitationally focused 
toward the NS. It is usually expressed as 

/?a = 2GMns/v^ = 3.7 X lO^^Vg-^ cm, (2.1) 

where vg is the wind velocity in units of 1000 km s^^ and we assumed that the orbital velocity of 
the star is negligible (Frank et al., 2002). Throughout the paper we fix the NS radius and mass at 
Rns=10^ cm and Mns=1-4 M©, respectively. The fraction Mcapt/M^ of the stellar wind mass loss 
rate (Mw) captured by the NS depends on Ra through (Frank et al., 2002) 

Mcapt/Mw ~ Rl/i^a^) = 2 X lO-^g-^fliod- (2-2) 

Here a=4.2x lO'^aiod cm is the orbital separation, aiod=PiodM3Q^, Piod is the binary orbital period 
in units of 10 days, and M30 is the total mass of the binary in units of 30 M© (we assumed circular 
orbits). 

- The magnetospheric radius, Rm, at which the pressure of the NS magnetic field {jx^l{%KR^^), 
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with n the NS magnetic moment) balances the ram pressure of the inflowing matter (Pw^^,)- In the 
case in which R]vi>Ra, the magnetospheric radius is given by (Davies & Pringle, 1981) 

Rm = 3.3 X lO^^MZ^^v^^^^a/^ln^l^ cm. (2.3) 

Here we assumed a non-magnetized spherically symmetric wind (Eisner & Lamb, 1977), with 
density^ p^{RM)^MJ{47ia^v^), a NS dipolar field with /i33=/i/10^^ G cm^, and M^(,=MJ\0^^ 
Mq yr^^ In the following sections we discuss the range of applicability of Eq. [2.3[ , and the regimes 
in which a different prescription for Rm should be used. 

- The corotation radius, Rco, at which the NS angular velocity equals the Keplerian angular ve- 
locity, i.e. 

7?co = 1 .7 X lO^^P^j^ cm. (2.4) 

Here Ps3 is the NS spin period in units of 10^ s. 

Changes in the relative position of these radii result into transitions across different regimes 
for the NS (Illarionov & Sunyaev, 1975; Stella et al., 1986). In particular, the accretion radius and 



magnetospheric radius depend on the wind parameters (see Eqs. |2JJ and [2.3| ), which can vary on a 
wide range of timescales (from hours to months). Therefore, variations in the wind parameters can 
cause the NS to undergo transitions across different regimes on comparably short timescales, thus 
opening the possibility to explain the properties of some classes of highly variable X-ray sources 
through them. Below we summarise the different regimes of a magnetic rotating NS, subject to a 
varying stellar wind. 

- Outside the accretion radius: the magnetic inhibition of accretion (RM>Ra). In sys- 
tems with RM>Ra the mass flow from the companion star interacts directly with the NS magneto- 
sphere without significant gravitational focusing, forming a bow shock at Rm (Harding et al., 1992; 
Toropina et al., 2006). At least in the front part of the shock, i.e. in the region around the stagnation 
point, the whole kinetic energy of the inflowing matter is converted into thermal energy. The power 
released in this region is of order 

Ashock ^ 2^mPwVw = 4.7 X 102^/?^ioV^£?io2m_6 erg s'^ (2.5) 

(R]yijo=RM/10^'' cm), and is mainly radiated in the X-ray band. We distinguish two different 
regimes of magnetic inhibition of accretion: 

• The superKeplerian magnetic inhibition regime: RM>Ra> Rco. In this case the magneto- 
spheric radius is larger than both the accretion and corotation radii. Matter that is shocked 
and halted close to Rm cannot proceed further inward, due to the rotational drag of the NS 
magnetosphere which is locally superKeplerian. Since magnetospheric rotation is also su- 
personic, the interaction between the NS magnetic field and matter at Rm results in rotational 
energy dissipation and thus, NS spin down. This process releases energy at a rate 

Lsd ^ nRliP^v^iRu^f ^ 3.7 x lO^^^ioM^e^rod^.s^ erg s"' , (2.6) 



which adds to the shock luminosity (Eq. 2.5). 



We approximated a-Rjyi— a , which is satisfied for a very wide range of parameters. 
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• The subKeplerian magnetic inhibition regime: Ra<RM<Rco- In this case the magneto- 
spheric drag is subKeplerian and matter can penetrate the NS magnetosphere through Kelvin- 
Helmholtz instabilities (KHI, Harding et al., 1992) and Bohm diffusion (Ikhsanov & Pustil'nik, 
1996). The mass inflow rate across Rm resulting from the KHI depends on the efficiency fac- 
tor t]kh~0.1, the shear velocity at Rm Vgh, and the densities pi and pe inside and outside the 
magnetospheric boundary at Rm, respectively. The luminosity released by accretion of this 
matter onto the NS is given by 

Lkh = 3.5 X lo3^T7KH/?Mio«rod^-6(Pi/Pe)'/'(l +Pi/Pe)"' erg s-\ {1.1) 

or 

Lkh = 8.8 x lo34T]KH/'s3'^Mio«ro>8"'^-6(pi/pe)'/'(l + Pi/pe)"' erg s-1, (2.8) 

if the shear velocity is taken to be the gas velocity in the post-shock region close to Rm or 
the rotational velocity of the NS magnetosphere, respectively (for simplicity, we consider 
throughout this paper Lkh equal to the largest of the above two values). The ratio Pi/pe can 
be calculated from the equation of the mass conservation across the KHI unstable layer, i.e. 

^MPeVconv - RuhPiVn{RM) , (2.9) 

where ht is the height of the unstable layer (Bumard et al., 1983). Throughout this paper we 
assume ht=/?M (Bozzo et al., 2008). The contribution of Bohm diffusion to the total mass 
inflow rate through the magnetosphere in the subKeplerian magnetic inhibition regime can 
be calculated according to Ikhsanov & Pustil'nik (1996), and we found that it is orders of 
magnitude smaller than that due to the KHI over the whole range of parameters relevant to 
this work. Similarly, the contribution to the total luminosity resulting from the shock at the 
magnetospheric boundary can be neglected in this regime. 

- Inside the accretion radius: RM<Ra- Once Rm is inside the accretion radius, matter flowing 
from the companion star is shocked adiabatically at Ra (this produces an almost negligible contri- 
bution to the total luminosity) and halted at the NS magnetosphere. In the region between Ra and 
Rm this matter redistributes itself into an approximately spherical configuration (resembling an 
"atmosphere"), whose shape and properties are determined by the interaction between matter and 
NS magnetic field at Rm (Davies et al, 1979; Davies & Pringle, 1981). A hydrostatic equilibrium 
ensues when radiative losses inside Ra are negligible^ ; the atmosphere is stationary on dynamical 
timescales, and a polytropic law of the form pocp^+^/" can be assumed for the pressure and den- 
sity of the atmosphere. The value of the polytropic index n depends on the conditions at the inner 
boundary of the atmosphere, and in particular on the rate at which energy is deposited there. Three 
different regimes can be distinguished: 

• The supersonic propeller regime: Rco<RM<Ra- In this case the rotational velocity of the NS 
magnetosphere at Rm is supersonic; the interaction with matter in the atmosphere leads to 
dissipation of some of the star's rotational energy and thus spin-down. Turbulent motions are 

^We checked this is verified for all the case of interest for tliis paper. 
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generated at Rm which convect this energy up through the atmosphere, until it is lost at its 
outer boundary. In this case n=l/2, and the balance between the magnetic and gas pressure 
gives 

Rm^23x lO'\%lMZfvf^f cm. (2. 10) 

Matter that is shocked at ~Ra, reaches the magnetospheric boundary at Rm where the inter- 
action with the NS magnetic field draws energy from NS rotation. According to Davies & 
Pringle (1981), this gives the largest contribution to the total luminosity in this regime 

Lsd = 2nRlip{RM)cl{RM) ^ 

5.4 X lO^'M^eawd^s'^Mioi^ + 16/?aio/(3/?Mio))'/' erg s-^ (2.11) 
In the above equation Rai()=10"^°Ra and Cs(RM)=Vff (Rm)=(2GMns/Rm)^''^- 

The subsonic propeller regime: RM<Ra> Rco> ^w<Mim- The break down of the supersonic 
propeller regime occurs when Rm<Rco> i-C-. when the magnetosphere rotation is no longer 
supersonic with respect to the surrounding material. The structure of the atmosphere changes 
and the transition to the subsonic propeller regime takes place. Since the rotation of the mag- 
netosphere is subsonic, the atmosphere is roughly adiabatic («=3/2), and the magnetospheric 
radius is approximated by (Davies & Pringle, 1981): 

Rm^2x \&^aJ^J4Zl'\^J^ pL^l'' cm. (2. 12) 

In the subsonic propeller regime, the centrifugal barrier does not operate because Rm<Rco, 
but the energy input at the base of the atmosphere (due to NS rotational energy dissipation) 
is still too high for matter to penetrate the magnetosphere at a rate Mcapt (Davies & Pringle, 
1981). Nevertheless a fraction of the matter inflow at Ra is expected to accrete onto the NS, 
due to the KHI and Bohm diffusion^. We found that, in all the cases of interest to this paper, 
the KHI provides the largest contribution to the total luminosity in the subsonic propeller 
regime. By using similar arguments to those above, we estimated in this case 

T -ISxIO^^n R' M (1 + 16W(5/?Mio))^/^(Pi/Pe)^/\ ^^_-i ,,13^ 

^s3afodV8(l+Pi/pe) 

The subsonic propeller regime applies until the critical accretion rate 

Miim_, = 2.8 X 102p^33a?odV8/?M^K,(l + 16/?aio/(5/?Mio))"'/^ (2.14) 

is reached, at which the gas radiative cooling (bremsstralhung) completely damps convective 
motions inside the atmosphere (Davies & Pringle, 1981). If this cooling takes place, direct 
accretion at a rate Mcapt onto the NS surface is permitted. 

The direct accretion regime: RM<Ra> Rco> Mw>Miim- If Rm<Rco and matter outside the 
magnetosphere cools efficiently, accretion onto the NS takes place at the same rate Mcapt (see 
Eq. 2.2) at which it flows towards the magnetosphere. The corresponding luminosity is 

Lace = GMNsMcapt//?NS = 2 X 10^^ M ^ea\^^v^^ erg s"^ ~ 2 x lO^^Mij erg (2.15) 

where Mi5=Mcapt/10^^ g s^^ This is the standard accretion regime; the system achieves the 
highest mass to luminosity conversion efficiency. 



To our knowledge this is the first apphcation of the KHI to the subsonic propeller regime. 



6 



Magnetars and SFXTs 



E. Bozzo 




Figure 2: (A) Upper panel: Variation of the luminosity through different regimes, as a function of the mass loss rate from the 
companion star. In this case the parameters of the model are fixed at ^^33=0.8, Ps3=l, and V8=1.2. Lower panel: Relative position of 
the magnetospheric radius, Rm (solid line), with respect to the accretion radius (dotted line), and the corotation radius R<;o (dashed 
line), as a function of the mass loss rate from the companion star. 

(B) Upper panel: Variation of the luminosity through different regimes, as a function of the mass loss rate from the companion star. In 
this case the parameters of the model are fixed at ^133 =0.001, Ps3=l, and vg=1.2. Lower panel: Relative position of the magnetospheric 
radius, Rm (solid line), with respect to the accretion radius Ra (dotted line), and corotation radius Rco (dashed line), as a function of the 
mass loss rate from the companion star. 



3. Transitions and paths across different regimes 

We explore here the conditions under which transitions across different regimes take place. 
As emphasised in § ^ these transitions occur when the relative positions of Rm, Ra, and Rco 
change in response to variations in the stellar wind parameters. In the following, since Rm depends 
only weakly on the orbital period and the total mass of the system, we fix aiod=l, and investigate 
variations in the other four parameters: JU33, Ps3, vg, and M_6- The equations that define the 
conditions for transitions between different regimes are 

Ojod magnetic barrier; (3.1) 
Ryi > Rco =^ Ps3 ^ 2.6M_y'*Vg ^^^a\l^iJ.]^^^ centrifugal barrier with Rm > R3,; (3.2) 

2/3 • —1/3 2/3 2/3 

Rm > Rco ^ Ps3 ^ l-^'^iod^-6 ^8 /^33 Centrifugal barrier with Rm < Ra', (3.3) 
■fs3 ^ 4.5M_g^^^^ajQj^Vg*'^^^/l33^^^ subsonic propeller ^ direct accretion. (3.4) 

As an example, in Fig. ^ we use the above equations to compute the luminosity swings for a sys- 
tem with jU33=0.8, Ps3=l, and V8=1.2. The lower panel of this figure shows that, for 0. 1<M_6<100, 
the magnetospheric radius crosses both the centrifugal (Rco) and magnetic (Ra) barriers. Corre- 
spondingly, the system moves from the superKeplerian magnetic inhibition regime, to the super- 
sonic and subsonic propeller regime, and, finally, to the direct accretion regime, giving rise to a 
six-decade luminosity swing from ~10^' to ~10^^ erg s^^ We note that a large part of this swing 
(about five decades) is attained across the transitions from the superKeplerian magnetic inhibition 
to the direct accretion regimes, which take a mere factor of ~5 variation of Mw 

Fig. ^ shows that, in the presence of a standard NS magnetic field (10^^ G), such abrupt lumi- 
nosity jumps are not expected for a very slowly rotating (1000 s) NS (the other system parameters 
are the same as those of Fig. |2|A), since the magnetospheric radius is smaller than both Ra and Rco, 
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Figure 3: (A) Upper panel: Variation of the luminosity through different regimes, as a function of the mass loss rate. In this case 
the parameters of the model are fixed at /l33=0.01, Ps3=0.1, and V8=1.2. Lower panel: Relative position of the magnetospheric radius, 
Rm (solid line), with respect to the accretion radius R^ (dotted line), and corotation radius R<;o (dashed line), as a function of the mass 
loss rate from the companion star. (B) Same as (A) but for ^133=0. 01, Ps3=0.1, and V8=2.2. 



for any reasonable value of M^. Therefore, the direct accretion regime apphes, with the luminosity 
proportional to M^. 

In Fig. ^ we show the transitions for a system with /i33=0.01 and Ps3=0.1. The wind velocity 
is V8=1.2 in Fig. ^A, and V8=2.2 in Fig. |3|B. These two figures show that, for sub-magnetar fields, 
a 100 s spinning NS can undergo a transition across the magnetic barrier (besides the centrifugal 
barrier), for suitable parameters (a high wind velocity in the case at hand). Such transitions take 
place over a more extended interval of mass loss rates. For instance Fig. ^ shows that an increase 
by a factor ~100 in the mass loss rate is required, in this case, to achieve a factor ^10^ luminosity 
swing; this is comparable with the magnetar case of Fig. 
Taking into account of the examples discussed above, we conclude that: 

- Long spin period systems (P^y^l) require magnetar-like B-fields (lJ-33^0.1) in order for a large 
luminosity swing (^10^) to arise from modest variations in the wind parameters (e.g. a factor ~5 
in M-f,)- These luminosity swings might result from transition across different regimes through both 
the centrifugal and magnetic barriers. 

- Shorter spin period systems (Ps^'^l) must posses lower magnetic fields (H3i,<^0.1) for similar 
transitions to take place. Somewhat larger variations in the wind parameters are required in order 
to achieve similar luminosity swings to those of the long period case, and transitions between dif- 
ferent regimes occur in most cases through the centrifugal barrier. 

- Few or no transitions are expected for systems with either high magnetic fields and short spin 
periods, or systems with lower magnetic fields and long spin periods. In the first case the centrifu- 
gal barrier halts the inflowing matter at Rm and accretion does not take place; such systems might 
thus be observable only at very low (X-ray) luminosity levels ('2:^10^^-10'^^ erg s^^). In the second 
case RM<Rcofor a wide range of wind parameters, accretion can take place, and a high persistent 
luminosity is released (ciilO^^-Kfi^ erg s^^). 



4. Application to SFXT sources 



In this section we propose that transitions across different regimes caused by relatively mild 
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variations of the wind parameters are responsible for the outbursts of SFXTs. As a case study 
we consider IGRJ17544-2619 (Sunyaev et al., 2003), a SFXT observed by Chandra during a 
complex transition to and from a ~1 hour-long outburst, yielding the first detailed characteriza- 
tion of a SFXT light curve over a wide range of luminosity. The spin period of IGRJ17544- 
2619 is presently unknown. In't Zand (2005) showed that four different stages, with very differ- 
ent luminosity levels, could be singled out during the Chandra observation: (a) a quiescent state 
with Lx— 2x10^^ erg s^^ (b) a rise stage with Lx— 1.5x10^^ erg s^^ (c) the outburst peak with 
Lx— 4x10^^ erg s^^, and (d) a post-outburst stage (or"tail") with Lx— 2x10^^ erg s^^ (see panel 
(a) of Fig. these luminosities are for a source distance of ~3.6 kpc, Rahoui et al., 2008). The 
maximum luminosity swing observed across these stages was a factor of >6.5 x 10^. 

Motivated by the evidence for >1000 s periodicities in XTE J1739-302 and IGR J 16479-45 14, 
we discuss first the possibiUty that IGR J17544-2619 contains a very slowly spinning NS. We use 
/i33=l, Ps3=1.3, V8=1.4, and show in Fig. ^(b) the different regimes experienced by such a NS as a 
function of the mass loss rate. ForM_6<20 the above values give RM>Ra and R]vi>Rco> such that 
superKeplerian magnetic inhibition of accretion applies. The expected luminosity in this regime, 
~10^' erg s^^ is likely outshined by the X-ray luminosity of the supergiant star (the companion 
star's luminosity is not shown in Fig. 0, but it is typically of order ~10^^ erg s^\ Cassinelli et 
al., 1981). We conclude that the lowest emission state (quiescence) of IGR J 17544-26 19 can be 
explained in this way, with the companion star dominating the high energy luminosity (In't Zand, 
2005). The rise stage is in good agreement with the subKeplerian magnetic inhibition regime, 
where the luminosity (~10^^ erg s^^) is dominated by accretion of matter onto the NS due to the 
KHI. The uncertainty in the value of h translates into an upper limit on the luminosity in this regime 
which is a factor of ~10 higher than that given above (Bozzo et al., 2008). During the outburst 
peak the direct accretion regime must apply at a mass loss rate of M_6=500. In this interpretation 
direct accretion must also be at work in the outburst tail at M_6~3, where a slight decrease in 
would cause the magnetic barrier to close and the source to return to quiescence. According to this 
interpretation, if IGR J17544-2619 has a spin period of > 1000 s, then it must host a magnetar. 

Panel (c) of Fig. ^ shows an alternative interpretation of the IGR J17544-2619 light curve, 
where we fixed /X33=0.08, Ps3=0.4, and vg^l. For this somewhat faster spin (and lower magnetic 
field), the luminosity variation is mainly driven by a transition across the centrifugal barrier (as 
opposed to the magnetic barrier). In this case, the quiescent state corresponds to the supersonic 
propeller regime (M-6 <0.6), the rise stage to the subsonic propeller (0.6<M_6<2), while both 
the peak of the outburst and tail take place in the direct accretion regime at M_6=200 and M_6=10, 
respectively. Assuming an even faster NS spin period for IGR J17544-2619, a weaker magnetic 
field would be required. In panel (d) of Fig. 0, we show the results obtained by adopting /X33=0.001 , 
Ps3=0.01, and V8=2. The '^lO^^ erg s^^ luminosity in the subsonic propeller regime compares well 
with the luminosity in the rise stage. However, the luminosity of the supersonic propeller regime is 
now significantly higher than the quiescence luminosity of ~10^^ erg s^' (this is a consequence of 
the higher value of for which the supersonic propeller regime is attained in this interpretation). 
We note that, the whole luminosity swing takes place for a wider range of mass loss rates, and the 
outburst peak luminosity requires M_6— 3000, an extremely high value even for an OB supergiant. 

Interpreting the properties of IGR J17544-2619 in terms of a NS with a spin periods <ClOO s 
is more difficult. For instance, for the subsonic propeller regime to set in, the mass loss rate 
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corresponding to the transition across Rm^Rco rnust be lower than the limit fixed by Eq. 2.14. 
If instead the transition takes place at higher mass loss rate, the system goes directly from the 
supersonic propeller to the direct accretion regime (or vice versa), bypassing the subsonic propeller: 
therefore, the rise stage would remain unexplained. Since fast rotating NSs require lower magnetic 
fields for direct accretion to take place while in outburst, Eq. 2.14 is satisfied only for very high 
wind velocities (v8>2-3). On the other hand, an increase by a factor of ~2 in the wind velocity 
(with respect to the longer spin period solutions) would give a substantially lower Mcapt, such that 
the subsonic and the direct accretion regime luminosities fall shortwards of the observed values 
(unless unrealistically high mass loss rate are considered). Based on the above discussion, we 
conclude that IGR Jl 7544-2619 likely hosts a slowly rotating NS, with spin period >100s. Whether 
the magnetic barrier or the centrifugal barrier sets in, causing inhibition of accretion away from 
the outbursts, will depend on whether the spin period is longer or shorter then ^1000 s. 

As another example we discuss the case of IGR J16465-4507, a SFXTs with a spin period of 
228 s. The luminosity behaviour of this source is still poorly known. An outburst at 5 x 10^^^ erg s^^ 
was observed with INTEGRAL (assuming a distance of 12.5 kpc, Lutovinov et al., 2005; Smith, 
2004), which did not detect the source before the outburst down to a level of 5x10^^ erg s^'. 
About a week later, XMM-Newton revealed the source at 5x 10^^ erg s^' and discovered the 228 s 
pulsations (Lutovinov et al., 2005). If the direct accretion regime applied all the way to the lowest 
luminosity level observed so far, then an upper limit of /i33~0.004 would be obtained by imposing 
that the NS did not enter the subsonic propeller regime. On the other hand, if the luminosity 
measured by XMM-Newton signalled that the source entered the subsonic propeller regime, while 
direct accretion occurred only during the outburst detected by INTEGRAL, then a considerably 
higher magnetic field of /X33~0.07 would be required. 

The above discussion emphasizes the importance of determining, through extended high sen- 
sitivity observations, the luminosity at which transitions between different source states occur, in 
particular the lowest luminosity level for which direct accretion is still at work. In combination 
with the NS spin, this can be used to infer the NS magnetic field. Alternatively accretion might take 
place unimpeded at all luminosity levels of SFXTs, a possibility which requires a very clumpy wind 
as envisaged in other scenarios (Walter & Zurita Heras, 2007). In this case the NS magnetic field 
can be considerably lower than discussed here. 



5. Conclusions 



In this paper we reviewed the theory of wind accretion in HMXBs hosting a magnetic NS 
with a supergiant companion, and considered in some detail the interaction processes between the 
inflowing plasma and the magnetosphere, that are expected to take place when direct accretion 
onto the NS surface is inhibited. We then applied this theory to SFXTs and showed that their large 
luminosity swings between quiescence and outburst (up to a factor of ~10^) can be attained in 
response to relatively modest variations of the wind parameters, provided the system undergoes 
transitions across different regimes. We found that such transitions can be driven mainly by either: 
(a) a centrifugal barrier mechanism, which halts direct accretion when the NS rotation becomes 
superKeplerian at the magnetospheric radius, or (b) a magnetic barrier mechanism, when the mag- 
netosphere extends beyond the accretion radius. Which mechanism and wind interaction regime 
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Figure 4: Application of the gated model to the IGRJ17544-2619 transition from quiescence to outburst. Panel (a): A schematic 
representation of the Chandra light curve of IGRJ17544-2619 obtained by using the segments that were not affected by pile-up (see 
both panels of figure 2 in In't Zand, 2005). The different luminosity stages are clearly visible. According to In't Zand (2005), the 
count rates on the y-axis correspond to 2 x lO'" erg s^', 1.5 x lO'^ erg s^', 2 x 10^^ erg s"', and 4 x 10^^ erg s^', in the quiescence 
state, the rise state, the tail, and the peak of the outburst, respectively. Panel (b): Interpretation of the quiescence to outburst transition 
of IGR J 17544-26 19 in terms of the magnetic barrier model. The dotted horizontal lines mark the luminosity that divide the different 
regimes. The parameters of the model are fixed atPs3=l-3, V8=1.4, and /i33=l. Panel (c): Interpretation of the quiescence to outburst 
transition of IGRJ17544-2619 based on the centrifugal barrier model. The parameters of the model are fixed at ^133=0.08, Ps3=0.4, and 
V8=l. Panel (d): Same as panel (c) but here the parameters of the model are fixed at /i33=0.001, Ps3=0.01, and V8=2. 



apply will depend sensitively on the NS spin period and magnetic field, besides the velocity and 
mass loss rate in the supergiant's wind. In particular, the magnetic barrier mechanism requires long 
spin periods (>1000 s) coupled with magnetar-like fields (>10^'* G). On the other hand, magnetar- 
like fields would also be required if the centrifugal barrier sets in at relatively high luminosities 
(> 10^^ erg s~^) in NSs with spin periods of hundreds seconds. 

Evidence has been found that the spin periods of a few SFXTs might be as long as 1000-2000 s. 
Motivated by this, we presented an interpretation of the activity of IGR J17544-2619 (whose spin 
period is unknown) in terms of the magnetic barrier by a 1300 s spinning NS and showed that 
the luminosity stages singled out in a Chandra observation of this source are well matched by 
the different regimes of wind-magnetosphere interaction expected in this case. We discussed also 
an interpretation of this source based on the centrifugal barrier and a sUghtly shorter spin period 
(400 s), which can reproduce the luminosity stages comparably well. We emphasise that in both 
solutions the required magnetic field strength (>10^^ G and >8xlO^^ G, respectively) are in the 
magnetar range. 

While the possibility that magnetars are hosted in binary system with supergiant compan- 
ions has been investigated by several authors (e.g., Zhang et al, 2004; Liu & Yan, 2006), clear 
observational evidence for such extremely high magnetic field NS in binary systems is still miss- 
ing. According to the present study, long spin period SFXTs might provide a new prospective for 
detecting and studying magnetars in binary systems. 
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